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Abstract 

We consider long strips of finite width L < 13 sites of ferromagnetic Ising 
spins with random couplings distributed according to the binary distribution: 
— " Jo) + ^{Jij ~ ''^Jo)), < r < 1. Spin-spin correlation func- 

tions < aoCR > along the "infinite" direction are computed by transfer-matrix 
methods, at the critical temperature of the corresponding two-dimensional 
system, and their probability distribution is investigated. We show that, 
although in-sample fluctuations do not die out as strip length is increased, 
averaged values converge satisfactorily. These latter are very close to the crit- 
ical correlation functions of the pure Ising model, in agreement with recent 
Monte-Carlo simulations. A scaling approach is formulated, which provides 
the essential aspects of the R- and L- dependence of the probability distribu- 
tion of In < CTQaR >, including the result that the appropriate scaling variable 
is R/L. Predictions from scaling theory are borne out by numerical data, 
which show the probability distribution of In < aoa^ > to be remarkably 
skewed at short distances, approaching a Gaussian only as R/L ^ 1 . 
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I. INTRODUCTION 



Most studies of random magnetic systems focus on whether or not quenched disorder 
destroys a sharp phase transition and, in the latter case, whether critical exponents are the 
same as for the corresponding pure magnets • Less attention has been paid to the un- 

derlying probability distribution functions which govern the behaviour of sample-averaged 
thermodinamic quantities, and which are expected to be universal in certain circumstances 
(see below). Early work on probability distributions of correlation functions concentrated, 
as numerical applications were concerned, on strictly one- dimensional systems The 
behaviour, under renormalization group transformations, of the distribution of e.g. con- 
ductivities in percolation-resistor networks 0, or interactions in spin glasses has been 
studied as well. More recently, the probability distributions of bulk quantities such as en- 
ergy, magnetisation, specific heat and susceptibility of disordered Ashkin- Teller models have 
been investigated in two dimensions by Monte-Carlo simulations. Bond distribution 
functions in one-dimensional quantum spin systems have been revisited very recently [pJ]] . 

Here we deal directly with spin-spin correlation functions on finite-width strips of two- 
dimensional disordered Ising systems. The basic motivation for using this geometry is the 



fact that strip calculations, in conjunction with finite-size scaling concepts are among 

the most accurate techniques to extract critical points and exponents for non-random low- 
dimensional systems [l^ll^ . The rate of decay of correlation functions determines correlation 
lengths along the strip. These latter are, in turn, a key piece of Nightingale's phenomenologi- 
cal renormalisation scheme |]T^|T5[, and have been given further relevance via the connection 



with critical exponents provided by conformal invariance concepts ||16[. Early extensions 



of strip scaling to random systems [T^ have since been pursued further 1 1^,1^, 2^ and put 



into a broader perspective. Though this has been done with the help of ideas arising from 
the study of probability distributions the behaviour of the probability distribu- 

tions themselves has not been closely investigated in strip geometries. In particular, their 
evolution toward the two-dimensional system's form as strip width increases has not been 
analysed. 

We consider a two-dimensional random-bond Ising model on a square lattice with a binary 
distribution of ferromagnetic interaction strengths, each occurring with equal probability: 

PiJ^J) = l{SiJ^J - Jo) + SiJ.j - r Jo)), < r < 1 . (1) 



For this case, the transition temperature Pc = ^/ksTc is exactly known from duality pT| , p2 



sinh(2/5cJo) smh{2(3^rJo) = 1 . (2) 

We have studied strips of width L < 13 sites, with periodic boundary conditions, and length 
= 10^ sites. Throughout this work we fix r = 1/4 and T = Tc (1/4) as given by Eq. 
(I). Numerically, Tc(l/4)/Jo = 1.239... (to be compared with Tc{l)/Jo = 2.269... ). 
Using this value of r ensures that disorder effects are rather strong, while at the same 
time one keeps a safe distance from the percolation regime at r = (near which crossover 
to geometry-dominated behaviour is expected to complicate the picture); this choice also 
coincides with that used in several recent Monte-Carlo simulations P,^, thus comparison 
(when appropriate) is made easier. The choice of T = Tc is important, as it is here that the 
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probability distributions are expected to have a non-trivial universal form; furthermore, the 
extensive literature on critical correlations for pure systems, both making explicit connection 
to conformal invariance ideas |]16| and previous to that (see e.g. Ref. [0] and references 
therein), is an important reference frame against which to set our results. 

In what follows, we first illustrate the role played by intrinsic fluctuations in the proba- 
bility distribution of correlation functions, and show that even though these do not die away 
for large samples, the sample-to-sample fluctuations of averaged values do go down as sam- 
ple size increases. Next we compare our results for averaged critical correlations with those 
for a pure system, in order to check on a recent proposal arising from Monte-Carlo data 



which implies equality, within error bars, of the corresponding quantities. We then go on to 
identify the key features of the shape of distributions, and investigate their variation with 
distance R and strip width L. A simplified scaling theory is formulated, which provides the 
essential aspects of the R and L- dependence. Numerical data for the probability distribu- 
tions of correlation functions bear out the main predictions of scaling theory, in particular 
the role played by the combination R/L as an appropriate scaling variable. 



II. INTRINSIC FLUCTUATIONS AND AVERAGES 

We calculate the spin-spin correlation function G (R) = < o"q(j)j >, between spins on the 
same row (say, row 1), and R columns apart, of strips with periodic boundary conditions 
along the vertical direction. This is done following the lines of Section 1.4 of Ref. |T^, with 



standard adaptations for an inhomogeneous system ||T9[. At each iteration of the transfer 
matrix from one column to the next, the respective vertical and horizontal bonds between 
first-neighbour spins are drawn from the bond probability distribution, Eq. (|l|). By shifting 
the origin along the strip and accumulating the corresponding results, one then obtains 
averages of the correlation function (or of any function F of it, such as its logarithm, which 
will be of particular importance in what follows), to be denoted by (G) (or {F (G)) ), the 
R dependence being implicitly understood. With strips of length = 10^ sites, we are 
able to produce 10^-10^ independent estimates of G for 7 < R < 100 which is the range of 
distances to concern us here. 

Normalized histograms of occurrence of the allowed values of the correlation function 
(or, rather more frequent below, of its logarithm) are produced by dividing a convenient 
interval of variation of G (or InG) into 10^ equal- width bins, and assigning each particular 
realization to the appropriate bin. For In G the interval ranging from In 10~^ to zero has 
proved generally adequate, except for R = 100 where the lower limit was pushed down to 
In 10"^^ 

For strictly one-dimensional disordered systems {i. e., chains) the average free energy 
(related to the largest Lyapunov exponent) has a normal distribution as befits a sum 

of random variables. Thus the fluctuations shrink with sample size (strip length) A^, and 
relative errors must die out as 1/ \fN . Correlation functions, on the other hand, are products 
of random variables, thus their distribution tends to a Zo(7-normal form as i? ^ cxd [^], that 
is, the probability distribution of InG approaches a Gaussian. However, the analysis of 
correlation functions turns out to be more complex than that of the free energy, even on 
chains; a primary reason for this is that while the latter quantity is self-averaging in the 
sense defined above, the former is not: the usual Brout argument EH] cannot be applied, as 
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explained e. g. in Ref. 0. The width of the probabihty distribution of correlation functions 
is then expected to be a permanent feature, which will not vanish (at least trivially) with 
increasing sample size. 

We have found that on finite-width strips the width of the distribution tends to stay 
essentially constant as varies. A graphic illustration is provided in Fig. |I] where the 
horizontal variable is In G, which turns out to be convenient for most purposes (see below) . 
Increasing simply smooths out the histogram; averages such as (G) or (InG) hardly 
move, the same being true of the width. This is easier to notice by comparing the Gaussians 
fitted to peak at (InG) and with width given by the root-mean-square deviation A(lnG) = 



{([InG-(lnG)]^)} 

Though neither A(lnG') or AG vanishes, it is still possible to extract valuable infor- 
mation from averaged values, the dispersion of which among independent samples (to be 
denoted respectively A(lnG) or A(G)) does shrink with increasing sample size. Fig. ^ shows 
the typical dependence of relative fluctuations, AG/{G) and A{G)/{G), with strip length 
A^. Varying the number n of distinct (that is, A^-lattice spacings long) samples within a 
reasonable interval, say n = 5-50, changes A(G) only slightly, consistent with a Xj ^fn- 
dependence to be inferred from standard arguments. From an investigation of A{G)/{G) 
for distances in the range R = 7-50 and strip widths up to L = 7, it turns out that both the 
order of magnitude and A^-dependence (roughly I/a/ZV) depicted in Fig. ^ are typical. The 
behaviour of A(lnG')/(lnG) is entirely similar. We can thus predict (see Fig. |^) that the 
fluctuations A(lnG)/(lnG) and A{G) / (G) will be of order 1% or just under that for strips 
of length A^ = 10^. This will be enough for our purposes here. Similar considerations have 
been used elsewhere in strip studies [Q, and seem to have been followed also in Monte-Carlo 
calculations of correlation functions in finite (L x L) systems [23|. 
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III. COMPARISON WITH PURE-SYSTEM CRITICAL CORRELATIONS 



It has been found in Monte-Carlo simulations |^3[ that the average correlation function 
at criticality of a random-bond Ising system is numerically very close to that for a pure 
system at its own critical point. Below we check on the corresponding quantities for the 
strip geometry. 

The spin-spin correlation function for the pure Ising model at T = Tc on a strip of width 



L is known from conformal invariance |16| to vary, for large R, L as: 



for spins along the same row as is the case here. The proportionality factor can be obtained 
from the exact square-lattice (L — oo) result p4|, < crla\ >= 0. 70338//?^/^. In Fig. ^ we 
show, for i? = 7 and 20, data for (G) and exp(lnG') together with a continuous curve for the 
pure system. The latter passes through numerically calculated points for L < 15 ( Eq. (|) is 
in error by one part in 10^ for L = 15, R = 7 and less than that for R = 20) and follows Eq. 
(0) for larger L. Using for the horizontal axis guarantees that the pure-system curve 
approaches the vertical axis linearly. However, it still shows high curvature even for the 
largest values of L attainable in our random-system calculations. This warns us to refrain 
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from extrapolating our data for L oo. Even so, we can learn from finite-width results 
that {G) behaves very closely to its pure-system counterpart. This is in line with previous 



findings |T9[ according to which in a random system the correlation length ^ to be used in 
the exponent-amplitude relation of conformal invariance, ^ = L/vrr] |T^, is that obtained 
from the average decay of (G) against R. Thus one gets a picture in which = 1/4 as for 
the pure system [|l9l, consistent with 7/z/ = 7/4 obtained e. g. from strip calculations of 



the average susceptibility for the random system and the scaling relation 7/z/ = 2 — 77. 

Of course the present result reaches further, as one could conceive of a scenario where 
the correlations would differ in the pure and random systems, but decay asymptotically 
as i? ^ 00 with the same rate (thus giving the same ^). In fact, the decay of exp(lnG') 
against R is not too dissimilar to that of {G) for moderate disorder, and for the finite 
values of L within reach of calculation. Only a systematic study of extrapolation trends 
as L — i> cx), covering different degrees of disorder, shows how the respective correlation 
lengths are essentially distinct [T^. The physical origin of this lies in that, on account of the 



properties of the probability distribution of G (to be seen in detail below), the most probable 
value exp(lnG') does not coincide with the average one, (G) 0,0]. Accordingly, it has been 



shown by field-theoretic arguments [26| and supported by numerical work |T9| that the most 



probable, or typical, correlation function decays as 

exp(lnG) ~i?-^/^(lni?)-^/« , (4) 

while the logarithmic corrections are washed away upon averaging for (G), resulting in a 
purely algebraic decay with rj = 1/4. 

Quantitative analysis of the results displayed in Fig. |^ shows that , considering the 
central estimates {G) the ratio Q = {G {R, L,r,Tc{r))) /G {R, L,l,Tc{l)) is, in all cases, 
within 1.01 — 1.03. With estimated error bars of order 1% as explained in the previous 
Section, a very small amount of overshooting seems to persist which does not follow a 
definite trend against R/L (see Fig. ^). Monte-Carlo data show the corresponding ratio 



approaching unity from below as lattice size L increases p3|, in the region R/L ^ 1. We 
cannot go far into that region, as the maximum strip widths within reach are not much 
larger than 10, and randomness effects are significantly distorted for small R ^ 5. 

Though we are not able at this point to advance a quantitative argument, it is plausible 
that finite-size effects manifest themselves differently in strip and square geometries; taken 



together with those of Ref. we interpret our data as evidence in favour of pure- and 



random-system critical correlations being in fact equal, at least for R, L ^ 1 and R/L <^ 1. 

In the next two Sections, we exploit the features of the probability distribution of G, and 
show that the variable R/L is indeed at least approximately the convenient one to describe 
several relevant aspects of the problem. 

IV. PROBABILITY DISTRIBUTIONS: SCALING THEORY 

A. Relevant parameters 

Our starting point is the result that, in one dimension and for large R the probability 
distribution of G must be log- normal ^,^^0. The same is expected to hold on strips 
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provided that R/L ^ 1. We seek for deviations from Gaussian behaviour as one moves 
away from this hmit. 

In Fig. 1^ we show normahzed histograms of occurrence of In G for fixed L = 5 and R = 7, 
20 and 50 . Though to zeroth-order one could say that all plots look similar to Gaussians, 
the semblance is reduced as R decreases. 

A quantitative measure of departure from Gaussian behaviour is the (dimensionless) 
skewness S, defined as [P^]: 

^^(i^^^y (5) 

for a distribution with mean (x) and dispersion a. Of course, for a finite number of real- 
izations of a given probability distribution S itself will be subject to fluctuations. In what 
follows we shall always quote S with two significant digits, which will allow us to discern 
trends while staying reasonably within reliable margins. For L = 5 and R = 7, 20, 50 and 
100 (the latter not shown in Fig. ^ one has respectively 5* = —0.67, —0.41, —0.24 and 
—0.19. We shall analyse the R— and L— dependence of S below; for the moment note that 
it approaches zero with increasing R, as expected, and is always negative. This is partly 
because of the many-channel, incipiently two-dimensional character of correlations on the 
strip: qualitatively, if there is at least one path of "strong" bonds between two spins, their 
correlation is significant, and is not much enhanced if there are more strong-bond paths, 
thus the peak at large In G with an abrupt cutoff above the maximum; on the other hand, 
configurations without any strong-bond paths at all are possible but with low probability, 
giving rise to the "tail" of very low In G values. This argument explains the general trend 
towards increasing \S\ for R/L < 1 as well (see below). An extreme example of this trend in 
shown in Fig. |^ for L = 13, R = 7 where S = —0.85. The negative skewness effect depends 
also on the variable against which histograms are plotted: here we use In G because our goal 
is to check on departures from a log-normal distribution, thus it is the skewness of this plot 
which matters in the context. For instance, a plot of the same distribution against tanh~^ G 
would have positive skewness. 

We choose to characterize the distribution of InG by three quantities, namely mean 
((InG)), dispersion, or root-mean-square deviation (A(lnG) ) and skewness {S). In other 
words, we assume that the probability distribution of correlation functions is satisfactorily 
described by perturbative corrections to a log-normal form. As shown below, this works well 
in the present case of ferromagnetic disorder. If frustration effects are present, such as e. g. 



in random-field p8[ or spin-glass systems, it may be necessary to take recourse to additional 
parameters, or even to adopt a different perspective. We shall not deal with this matter in 
the present work. 



B. Scaling Theory 



The aspects just described are consistent with a scaling description given in this subsec- 
tion. It leads to further specific predictions regarding the form of distribution functions and 
their dependences on the variables R and L, which will be discussed subsequently in the 
light of the numerical results. 
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The approach combines two main features: (i) the appearance at large scales L, R ^ 

1 of universal aspects related to fixed point Hamiltonian and probability distributions of 
the disordered two-dimensional {2d) Ising model; and (ii) the crossover for R > L to 
width-limited behaviour characteristic of the one- dimensional (Id) version of the large scale 
universal properties. 

The procedure is to first scale n times by (length) rescaling factor b, where b"' = L, which 
takes the system to an equivalent linear chain. This step involves the scaling of joint probabil- 
ity distributions for the appropriate variables. If one starts from the critical condition of the 
random 2d Ising system irrelevant variables scale away, and one approaches asymptotically 
the fixed point Hamiltonian and distribution describing the universality class containing the 
2d random Ising system. This involves a minimal set of relevant random variablesjtj-"'' } (af- 
ter n scalings) and their universal probability distributions. The correlation function scales 
as follows: 

GL{R,{U})=b-''GL/biR/b,{t^^}) = ■■■ = L-''G^iR/L,{t^}) {n = \nL/\nb) . (6) 

We now have an equivalent Id system with Hamiltonian close to the fixed point Hamiltonian 
of the disordered 2d system. Since the correlations in a Id system are transmitted through 
each intermediate space point, a factorisation of Gi {R/L, {t'f^}) is suggested li R/ L > 1. 
This factorisation is into R/L factors corresponding to successive Lx L blocks of the original 
system [i.e. to single bonds of the renormalised Id system), labelled by s = 1, . . ., R/L; 
then Eq. (|) becomes 

R/L 

GUR,m=L-' nG'i(i,{ti"^K) . (7) 

s=l 

Here, {t\"'^}s are the renormalised random variables for the block labelled by s. It will 
be assumed later that these are largely uncorrelated from one block s to another, since 
the blocks were initially non-overlapping. This, and the other approximations leading to the 
approximate form Eq. (|^) will be tested by later comparison to the numerical results. Eq. (^ 
generalises a result for the pure case. There, for R/ L large, Gi (R/ L, . . .t* . . .) oc (A2/ Ai)^/^, 
where A2/A1 is a universal ratio of eigenvalues of the transfer matrix of the universal fixed 
point Hamiltonian of the pure 2d Ising class. Then 

GLiR,t*)ocL-^ exp[-R/^L] , (8) 

where = -^ln(Ai/A2) (= Tcq/L |[T^). It is perhaps interesting that the simplest b = 

2 Migdal-Kadanoff real-space renormalization group transformation gives Gi {R/ L,t*) = 
(t*)«/^, t* = 0.544, hence it gives 1] ~ 0.19. This unsatisfactory representation of the 
universal 77 in terms of a non-universal t* is due to not having allowed the Hamiltonian to 
adopt its universal form. 

We now return to the disordered case and consider the development of the probability 
distributions during the rescalings leading to Eq. (^. To allow for correlations it is necessary 
to consider the probability distribution Pt{ti} for the whole set of tj's. This is labelled by 
a parameter t setting the scale for all the tj's. The scaling of the distribution is given by a 
mapping Wb'. 
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P^;^ {U} = W, {P« {U}} = Wl [Pf,, {U}] , (9) 

where /, Z + 1 label two successive steps, and Wl denotes / iterations of the map. The 
parameter t also scales according to a renormalization group transformation characteristic of 
the 2d random Ising model. At the fixed point t* of that transformation, after n = InL/ ln6 
scalings with n large, will have produced a distribution P[^^ close to the universal 
invariant distribution P^t of the random 2d Ising model, which satisfies the fixed point 
equation 

p:,{t?i = w,{Pi,{u}} . (10) 

Then, employing Eq. (|^) and taking logarithms to obtain a sum of random variables on the 
right-hand side we find that at t*, after n scalings, the probability distribution P(q;) for 
InG^ {i.e., for the probability that InG^, {R, {tj}) takes the value a) is given by: 

I da exp{f3a)V{a) = I [Hdu] p/.") {U} exp{-f3r] In L) [] exp(/51nGi(l, {U},) . (11) 

s=l 

If n is large, as well as P/"'' ~ P/, , we have Gi{l, {ti}s) close to a universal function charac- 
teristic of the random 2d Ising class, since the renormalised bond regions labelled by each s 
are each a composition of many original bonds, at t*. 

To the extent that the probability distribution P/"^ factorises into parts q^*^ {ti}s corre- 
sponding to the different s's (to be tested) the result Eq. (p!T|) reduces to a form corresponding 
to the probability of a sum of random variables (InGi): 

J da exp{(3a)V{a) = exp{-(3r]\nL) (/(")(/5))^^'^ , (12) 

with 

I^'^KP) = I {Udu)^ g;r){t,},exp(/51nGi(l,{t,},)) , (13) 

where the subscript s indicates that all tj's are in the region corresponding to a given s. 

We now explore the consequences of the above results, and in particular the expected uni- 
versality of ql* , Gi{l, {ti}s) for L (= 6") large, for the distribution V{a) for InGiiR, {ti}s) ■ 

Eq. (|l^) shows that V {a — rj In L) corresponds to the probability distribution for a sum of 
R/L independent random variables (In 6*1(1, {ti}s)), and has the consequence that, if R/L 
is large, V{a) approaches a Gaussian distribution with mean, width and skewness given by: 

(a) = (^^^ m - r]\nL (14) 
{{a-{a))Y'=(^y^\ (15) 



L 



((a-(«))3)/((a-(a))y/2=(^)"'^%. (16) 
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The quantities m, w, s are characteristics of /*^"^(/?), which is related as follows to the 
(non-Gaussian) distribution function Q^^\a) for the logarithm of the nearest-neighbour 
correlation function of the (n-times rescaled) disordered 2d Ising system: 

/(")(/5) = Jda exp(/5a) Q!"^(a) , (17) 



with 



Qt\^) = J {lidti)^ qi?^{U}J{a-\nG^il,{U}s)) . (18) 

For L = large, all these quantities and therefore m, w, s will become universal (charac- 
teristic of the 2d random Ising class). 

The universal character at large L of the distribution V{a) for InG^, (R) at t*, its Gaus- 
sian form at large R/L and the results (|T4|)-(p!6|) and their interpretation above are the main 
conclusions of this Section. Note (a) that these conclusions include the result of Ref. ||^ for 
the Id case and (6) that the analogue of m for the pure case is the universal constant —ttt) 
(while w, s are zero). 

The essential points of the general discussion given above can be explicitly illustrated in 
the simple scenario provided by the renormalization group approach (blocking/decimation), 
allowing for just the random variables {tj = tanh/5Jj}. This forces the Hamiltonian to 
remain of Ising form (so the reservation expressed under Eq. (^ applies). We allow for spin 
rescaling, and for the distribution function g{ti) for each to develop towards fixed point 
universal form under the scaling [§,0, but we ignore correlations. Then Gi(l, {tj}) is just 
the renormalized variable t\^\ The scaling of ti will be of the form 

t^^) = R,{t?} (19) 

where the right-hand side is a function of Nf, variables tf ■* comprising the block. Then the 
distribution g^''\tf^) scales according to the following simplified version of Eq. (P): 




and t'''-' t^'"*"^) is the resulting change of scale of the distribution g. No scale change occurs 
if the initial distribution is set at the critical value of t^^\ Then, for large I, Eq. (^0]) gives 
the asymptotic approach to the universal fixed point distribution, from which m, w, s can 
be obtained. For an adequate description of this sort, the transformation Eq. ([19|) should 
give the proper zero value of the exponent a for the pure case: via the Harris criterion this 
makes the disorder marginally relevant and ensures that the width w does not scale away. It 
requires A*";, = A^, where A;, = b^^'^ is the eigenvalue of the pure version of Eq. (0), linearized 
about its fixed point. Procedures of this sort give (i) m not very different from its pure 
value; (m) negative skewness s. 



V. NUMERICAL RESULTS AND CONTACT WITH THE SCALING THEORY 

We begin by recalling that the results exhibited in Fig. |^ point out the importance of 
intrinsic widths in the critical random system {i.e. at t*). This is a central feature of the 
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scaling theory, through the appearance of universal distributions. Secondly, we recall that 
Fig. ^ provides evidence for Gaussian distributions for InG at large R, with narrowing 
relative widths as R increases. This is again a prediction of the scaling theory (see Eq. (|T^ 
and the discussion preceding it). 

The results already presented in Fig. |^ show that the overall dependence of exp(lnG') on 
L for fixed R approximately mimics that of (G), apart from a proportionality factor. The 
latter is, in turn, numerically very close (see Fig. ^) to that of the pure system, given to 
good approximation by Eq. (^. The proportionality factor is, however, R— dependent, as 
illustrated in Fig. this illustrates that, though both quantities decay exponentially with 
R as befits an essentially one-dimensional system, their respective correlation lengths differ, 
with well-known consequences [|^,[7i[T9| . The hnear dependence of (InG) = a against R is 
consistent with the scaling result (|14|). 

Also, at large R, the slope of the measured a versus R line is proportional to l/L, in 
agreement with the {R/L)m term dominant in Eq. (|T^); and the coefficient is consistent 
with having m not far from its pure value. The comparison of numerical results for (In G) 
and (G) with corresponding pure forms is given in Fig. |^ where Gpure is given by Eq. (^ 
above. 

A further, rather stringent test of the scaling predictions is the plot in Fig. ^of numerical 
data for A(lnG) against R/L. It can be seen that the (R/LY^'^ dependence Eq. (|15|) is 
reasonably followed in the large R/L regime where it was derived; and the data collapse of 
results for different R, L suggests that R/L is (as predicted) the appropriate scaling variable 
in the regime of large R, L independent of their ratio. The curve also gives evidence of the 
crossover to the universal width of the 2d random Ising system for L ^ R^ 1. 

In order to further test the scaling theory, and the suggestion that R/L is the apppro- 
priate scaling variable, we show in Fig. ^ numerical results for the skewness against R/L. 
Again the data collapse is satisfactory. The prediction of (R/L)^^^"^ behaviour (Eq. ([T6|)) 
for large R/L is again observed, and the crossover towards universal 2d behaviour is seen 
for L > i? > 1. 

VI. CONCLUSIONS 

We have studied properties of the probability distributions of correlation functions on 
finite-width strips of the two-dimensional random-bond Ising model at criticality. We have 
shown that even though intrinsic fluctuations in the probability distribution do not die away 
for large samples, the sample-to-sample fluctuations of averaged values do go down approx- 
imately with the square root of sample size as the latter increases. Results thus obtained 
for averaged critical correlations have been compared with those for a pure system, and we 
have found that the values of averaged correlations (G) are very close to the correspond- 
ing pure-system ones, consistent with recent Monte-Carlo data [^. The key features of 
the shape of distributions have been identifled, and a simplifled scaling theory has been 
formulated, which provides the essential aspects of the R and L- dependence. Numerical 
data for the probability distributions of correlation functions bear out the main predictions 
of scaling theory, in particular the role played by the combination R/L as an appropriate 
scaling variable. 
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We expect the approach outhned above, which consists in describing the probabihty 
distribution of InG by perturbative corrections to a log-normal form (thus characterized 
by only three quantities, namely mean ((InG)), width (A(lnG) ) and skewness (S)) to be 
appropriate in the present case of ferromagnetic disorder; it remains to be checked whether 
additional parameters, or even a change of perspective, will be necessary if frustration effects 
are present, e. g. in random-field pSf or spin-glass systems. We plan to undertake this task 
ClS db continuation of the present work. 
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FIGURES 

FIG. 1. Normalized histograms of occurrence of InG for L = 5, R = 20. (a) : N = 10^; (6) 
: AT = 10*^. Full vertical arrows at (InG); broken vertical arrows at ln(G) . Curves are Gaussians 
fitted to mean and root-mean-square deviation of InG, as calculated from respective realizations. 
Here, and in all subsequent figures, r = 1/4 and T = Tc{r). 

FIG. 2. Relative fluctuations within sample, AG/{G), and between sample-averaged values, 
A{G) / (G) , against strip length N. L = 7, R = 20, number of samples n = 20. 

FIG. 3. Averaged correlation functions for L = 5, 7, 9, 11, 13. (a) : R = 7; (6) : R = 20. 
Triangles: (G); squares: exp(lnG). Errors as defined in Section II. Continuous line: correlation 



functions for pure Ising model on strips |1£524| . 



FIG. 4. Ratio Q = {G {R, L,r,Tc{r))) /G {R, L,l,Tc{l)) against R/L for points of Fig. |. Error 
bars represent estimated fluctuations of order 1%. 

FIG. 5. Normalized histograms of occurrence of InG for fixed L = 5 and i? = 7, 20 and 50 
. Full vertical arrows at (InG); broken vertical arrows at ln(G) . Curves are Gaussians fitted to 
mean and root-mean-square deviation of In G, as calculated from respective realizations. 

FIG. 6. Normalized histogram of occurrence of InG for L = 13 and R = 7. Full vertical arrow 
at (InG); broken vertical arrow at ln(G) . Curve is Gaussian fitted to mean and root-mean-square 
deviation of InG, as calculated from respective realization. Skewness = —0.85 . 

FIG. 7. Semi-log plot of decay of (InG) and ln(G) against distance for fixed L = 5. Broken 
line gives slope as predicted by conformal invariance (L/^ = vrr/) with rj = 1/4. 

FIG. 8. Log-log plot of ratio between averaged correlation functions and Gpure against R/L, 
the latter as given in Fig. ^ [16,^. Triangles: data from (G); squares: data from (InG). 



FIG. 9. Log-log plot of width A(lnG) against distance R/L . Line has slope 1/2 . 
FIG. 10. Log-log plot of negative skewness —S against R/L . Line has slope —1/2 
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